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NYO-9^92 
Abstract 

Adlam  and  Allen  [1]  and  Davis,  Lttst,  and  SchlUter  [k] 
have  studied  nonlinear  plane-waves,  propagating  normal  to 
the  magnetic  field,  in  a  cold  plasma.   One  solution  of 
particular  Interest  is  a  solitary  wave,  or  single  pulse.   We 
present  a  method  for  solving  the  analogous  problem  for  a 
plasma  with  finite  temperature,  in  the  limiting  case  where 
the  amplitude  of  the  wave  is  small  --  and  where,  consequently, 
the  v/idth  of  the  wave  is  very  large.   Our  result  is  that  in 
a  v;arm  plasma  the  solitary  wave  of  small  amplitude  has  the 
same  form  as  in  a  cold  plasma,  but  has  a  width  which  is 
larger  than  its  width  in  a  cold  plasma  by  a  factor  which  is 
approximately 


Here  p  is  the  ratio  of  material  pressure  to  magnetic  pressure 
in  the  undisturbed  plasma;  and  m  ,  m_  are  the  respective  masses 
of  a  positive  ion  and  of  an  electron.   Evidently  this  factor 
is  appreciable,  unless  p  is  very  small. 


-  2 


NYO-9^92 


Table  of  Contents 


Page 


Abstract  2 
Section 

1.  Introduction  4 

2.  Equations  of  Motion.   Dlmenslonless 
Variables  and  Scaling  Transformations  6 

3-   Series  Solution  of  the  Collision-free 

Boltzmann  Equations.  13 

4.   Computation  of  the  Moments  19 

5-   Differential  Equation  for  B  22 

6.   Width  of  the  Solitary  Wave  25 

References  26 


-  5 


NYO-9^92 

The  Effect  of  Temperature  on  the  Width  of  a  Small -Amplitude, 
Solitary  Wave  in  a  Collision-Free  Plasma 

1.  Introduction. 

Adlam  and  Allen  [1]  and  Davis,  LOst,  and  Schlfliter  [4]  have 
studied  the  nonlinear  steady  plane  waves  which  can  propagate 
normal  to  the  magnetic  field  in  a  collision-free  plasma.  These 
authors  assumed  a  cold  plasma  --  in  which  the  temperature  of  the 
positive-ion  gas  is  zero,  and  the  temperature  of  the  electron 
gas  is  zero.  Thus  at  any  point,  in  a  cold  plasma,  all  ions  have 
the  same  velocity,  and  all  electrons  have  the  same  velocity.   A 
one -dimensional  steady  flow  of  a  cold  plasma  is  therefore 
governed  by  ordinary  differential  equation.  These  equations 
were  solved  by  Adlam  and  Allen,  and,  in  fuller  detail,  by  Davis, 
Liist,  and  Schlflter.   In  general  the  solution  represents  a  peri- 
odic wave  train  (which  is  not  sinusoidal,  in  this  nonlinear 
theory).  A  particular  solution,  of  special  Interest,  represents 
a  single  pulse,  or  solitary  wave.   Gardner,  et  al.  [5],  Rose  [14], 
Banos  and  Vernon  [2]  and  Vernon  [15]  have  obtained  similar 
solutions  from  various  two-fluid  theories  (based  on  suggestions 
due  to  Grad,  [7]  and  [8])  which  include  pressure  and  temperature. 

Our  object  is  to  investigate  the  effect  of  finite  plasma 
temperature  on  the  solitary  wave  solution.   If  no  approximations 
are  made,  (other  than  that  the  plasma  is  collision-free),  the 
equations  of  motion  are  too  difficult  to  be  solved  for  a  warm 
plasma  (with  temperatiore ) .  We  are  able,  however,  to  find  a  solu- 
tion in  the  limiting  case  of  a  wave  of  small  afnplitude.   Our 
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approximation  here  does  not  lead  to  a  linear  theory;  a  linear 
theory  could  never  yield  a  solitary  wave  as  solution.   The 
reason  is  that  the  width  of  the  solitary  wave  tends  to  infinity 
as  the  amplitude  tends  to  zero;  so  that  the  Independent  variable 
must  be  scaled  accordingly.   In  fact,  we  assume  the  same  scaling 
laws  for  the  warm  plasma  as  obtains  for  the  solitary  wave  in  the 
cold  plasma.   That  is,  vie   assume  that  the  width  of  the  wave 
varies  inversely  as  the  square  root  of  its  amplitude.   This  was 
shown  for  the  cold  plasma  by  Gardner  and  Morikawa  [6].   Thus 
a  length  over  which  the  magnetic  and  electrjc  fields  vary 
appreciably  is  very  large  in  this  limit.   Consequently,  one  can 
obtain  a  formal  solution  of  the  collision-free  Boltzmann  equation 
by  a  successive-approximation  calculation.   This  calculation  yields 
expressions  for  the  distribution  functions  of  the  ions  and  elec- 
trons in  terms  of  the  electric  and  magnetic  field  components  and 
their  derivatives.   Inserting  these  expressions  into  Maxwell's 
equations,  and  retaining  only  the  terms  of  lov;est  significant 
order,  v;e  obtain  a  differential  equation  for  the  asymptotic  shape 
of  the  v/ave.   This  equation  has  the  same  form  as  the  analogous 
equation  for  the  cold  plasma;  however  the  coefficients  depend  on 
the  plasma  temperature.   The  final  conclusion  which  can  be  drawn 
is  that  the  effect  of  temperature  is  to  multiply  the  width  of  the 

_  . 

We  point  out  that  the  inherent  non-linearity  (even  for  small 
amplitude)  of  the  solitary  wave  solution  in  the  study  of  gravity 
waves  in  water  of  finite  depth  has  been  knov/n  for  many  years 
(cf.  Korteweg  and  de  Vries,  [10],  and  Keller,  [9])«    Other 
physical  problems  which  exhibit  this  peculiarity  are  described 
more  recently  by  Morlkav-a  [12],  [13]' 


Adlam-Allen  solitary  wave  by  a  factor  which  is  approximately 

(assuming  for  simplicity  that  the  ion  temperature  and  electron 
temperature  are  equal).   Here  p  is  the  ratio  of  material  pres- 
sure to  magnetic  pressure  in  the  undisturbed  plasma;  and  m  is 
the  mass  of  a  positive  ion,  and  m  the  mass  of  an  electron. 
Evidently  this  factor  is  appreciable  unless  p  is  very  small. 
Another  way  of  stating  the  result  is  that  the  cold-plasma  treat- 
ment of  the  solitary  wave  is  correct  only  if  the  electron  sound  speed 
is  small  compared  to  the  Alfven  speed  in  the  undisturbed  plasma. 
This  result  may  be  pertinent  to  the  study  of  weak,  collision-free 
plasma  shocks  (cf.  Morawetz,  [11]). 


2.   Equations  of  Motion.   Dimensionless  Variables  and  Scaling 
Transformations 
Our  problem  is  to  study  a  plane  wave,  propagating  with  speed 
U^  in  the  negative  x  direction,  in  a  plasma,  initially  in  equilib- 
rium, and  permeated  by  a  constant  uniform  field  of  strength  B  , 
in  the  z   direction.   We  shall  formulate  this  problem  in  a 
coordinate  system  which  moves  with  the  wave.   We  assime  that  in 
this  coordinate  system  we  have  a  steady-state  flow  of  the  plasma. 
Due  to  the  coordinate  transformation  there  is  now  a  constant 
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electric  field,  in  the  y  direction,  of  magnitude  E  =  B  U  . 
We  use  m.k.s.  electrical  units.   Here,  of  course,  B  means  the 
undisturbed  magnetic  field  as  it  appears  in  the  moving 
coordinate  system.   The  plasma  consists  of  positive  ions  and 
electrons  continuously  distributed.   We  assume  that  the  motion 
of  an  ion  or  electron  is  determined  only  by  the  force  exerted  on 
it  by  the  macroscopic  electric  and  magnetic  fields.   The  magnetic 
field  is  entirely  in  the  z"  direction,  and  has  strength  B(x) 
(depending  only  on  x,  and  independent  of  time,  in  the  frame  of 
reference  fixed  to  the  wave).   The  electric  field  has  the  com- 
ponent E  =  B  U   and  also  the  component  E  (x ) ,  due  to  charge 
separation.   The  plasma  proper  is  described  by  the  distribution 
functions  of  the  ions  and  electrons,  which  are  time-independent, 
and  depend  on  only  the  one  space  coordinate  x. 

Let  u,  V  denote  the  x  and  y  components  of  the  velocity  of 
a  particle.   Let  T  (x,u,v)  and  r_(x,u,v)  denote  the  density  of 
the  ions  and  the  density  of  the  electrons  in  the  x  u  v  phase 
space.   Then  the  distribution  functions  f  ,  f"_   satisfy  the 
following  collisionless  Boltzmann  equations: 

(1)       ^-Z^  +  l-   (E^  +  Bv)  -±  +  I-  (E  -  JBu)  ^  =  0 


(2)       u  — :i  -  ^  (E  +  Bv)  — ^  -  |-  (E  -  Bu)  -3^  =  0 
^x    ^-       ^  ^u    "^-   ^      ^v 


Here  e  is  the  electronic  charge,  and  m  is  the  mass  of  an  ion. 
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and  m_  the  mass  of  an  electron.   The  electric  and  magnetic 
fields  are  governed  by  the  Maxwell  equations. 


^2^     o      _  -  _  _ 

(3)       -i^  =  M-C^e  [//  f  dudv  -  //  f  dudv] 
dx  "^  " 


(2|]       ^  =  _p.e  [//  vf  dudv  -  //  vf  dudv] 
dx  "^  " 


where  [.i   is  the  permeability  of  space  and  c  the  speed  of  light. 
(We  use  M.K.S.  units.) 

We  are  interested  in  solutions  which  tend  to  a  constant 
state  as  x  tends  to  -oo  .   That  is,  we  suppose  the  following 
conditions : 

For  X  — ^  -00  , 


(5)       f^(x,i5,v)  ->  -^^((u-U^)^  +  v^) 


(6)       f_(x,u,v)  ->  ^_((u-U  )2  +  v^) 


(7)  B(x)  -^   Bq   (a  constant) 

(8)  E^{x)  -^   0 

where 

(9)  E  =  B  U  , 
V  ^  /        y    o  o ' 

and  where  the  number  densities  of  ions  and  electrons  are  equal 
//  T^  dudv  =  n  =  //  ^_dudv 
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Here  W    and  W    are  the  undisturbed  distribution  functions  for 
ions  and  for  electrons.   The  most  interesting  case,  of  course, 
is  the  case  wherein  ?",  ,i/     are  Maxwellian  with  equal  temperature 
(not  a  necessary  restriction). 

It  is  easily  verified  that  one  solution  is  given  by 


+  =  ^+>  7_  =  f_,   B(x)  =  Bq,   E^  = 


This  is  the  trivial  solution  --  corresponding  to  no  disturbance 
of  the  plasma.   It  is  not  really  clear  than  any  other  steady 
state  solution  does  in  fact  exist.   However,  we  shall  assume 
that  such  a  non-trivial  solution  exists,  and  furthermore  that 
it  depends  asymptotically  on  a  parameter  €  in  a  certain  way.   We 
then  proceed  to  find  the  solution  for  small  e,  correct  to  lowest 
order  in  e . 

To  define  e,  we  first  must  define  a  few  other  parameters 
of  importance.   Let  p  ,  p  denote  the  initial  pressures  of  the 
ion  gas  and  of  the  electron  gas;  then 


(10)  P+  =  I  m^  //[(u-Uq)^+  v2]^^dudv 

(11)  P_  =  I  ni_  //[(u-Uq)^+  v^]^_dudv   . 

Let  a  be  the  sound  speed,  and  A  the  Alfven  speed,  of  the 
undisturbed  plasma;   then 

,   ,  2   2[p  +  pj 

(12)  a^  =  r^ , 

^Qlm^  +  m_] 
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(13)  A^  = 2 

M-nQ[m_|_  +  m_] 

Let   the  nuraber  M  be   defined  by 

(14)  M^  =  ° 


""5 2 


Then  M  is  akin  to  Mach  number,  being  the  ratio  of  the  flow 
speed  Uq  to  the  speed  (A  +  a"")    ,   v/hlch  is  the  propagation 
speed  of  a  small-amplitude  hydroraagnetic  wave  normal  to  the 
magnetic  field.   We  define  e  by  setting 


(15)  M^  =  1  +  €^. 


We  define  dlmensionless  variables  (without  bars)  as 
follows : 

(16)  x=i     +Ox 

^      '  e      eB^ 

o 

(17)  u   =   Uq(1    +  u)     ;      V   =  U^v 

(18)  B  =  B^(l   +   e^B) 

(19)  E^  =  e\u^E 

(20)  r   =  ^  f    .    r   .^  %  f 

n^  ^  ^  n 


(21)  f     =  ^f   (u-^  +  v^)}     W     =     o^(u^  +  v'^) 

The  reason  for  introducing  the  powers   of  e   into    (l6),    (l8). 
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(19)  is  that  we  assume  that  the  length  scale  of  the  disturbance 

is  of  order  lA,  and  that  the  disturbance  in  B  is  of  order  e  , 

3 
and  the  disturbance  E  is  of  order  e    ,   as  e  tends  to  zero.  These 

scaling  rules,  which,  as  we  shall  see,  are  self-consistent,  were 

motivated  by  their  validity  for  the  solitary  wave  in  a  cold 

plasma.  The  explicit  solution  is  available  for  the  solitary 

wave  in  a  cold  plasma.   And,  in  agreement  with  the  above,  it  is 

found  that  the  non-trivial  solution  exists  only  if  M  >  1;  and 

for  M  =  1  +  e  ,  and  E  is  of  order  e-^,  (as  shown  by  Gardner  and 

Morikawa ) . 

In  terms  of  the  new  variables  given  by  (l6)-(20),  the 

Boltzmann  equations  take  on  the  following  forms : 

Sf.    p  Sf.         Sf.    2  Sf.    ,  Sf. 
(22)      ^  +  c^B  ^  =  €(l+u)  ^  -  e^B  ^  +  c^E  ^ 

Sf     p  Bf  Sf     P  Sf     ,  Sf 

(25)      ^  +  e  B  ^  =  -ea(l+u)  ^  -  c  B  ^  -f  e^E  ^  . 


Here  0  is  a  polar  coordinate  in  the  u,v  plane;  v/e  define  w,  6   by 


(24)  w  =  u^  +  v^ 

(25)  u  =  w-'-^^cos  e  ,   V  =  w-^^^sin  0 
so  that 

The  coefficient  a  in  (  25)  is  the  ratio  of  the  masses  of  the 
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tv;o  types  of  particle: 
(27)      a  =  jjj^. 

The  Maxwell  equations  now  take  on  the  following  forms 


hf^  ,  V  A^  dE 


(28)  €"(l+a)  ^  SI  =  [//  f  dudv  -  //  f  dudv] 

(29)  e^[-(l+a)+2(l+€^)(P^+aP_)]g  =  (l+e^)  [//vf_^dudv 

-  //  vf_dudv]. 

Here  P  and  P_  are  dlmensionless  pressures  defined  by 

(50)      P  =  ±-^     ,        P   =  ^   . 

n  m  n""  n  m  U 

o  +  o  o  -  o 

Finally,  we  have  the  following  conditions  as  x  tends  to 
—CO  ; 

For  X  — >  -00, 

(31)  f^(x,u,v)  -^^'^(u^  +  v^) 

(32)  f_(x,u,v)  -^  •^_(u^  +  v^) 

(33)  B(x)  -»0 

(3^)      E(x)  ->  0 

and 

(35)      //  "^^dudv  =  1,    //  V^_dudv  =  1 
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(36) 


v_ 


I  //  (u^  +  v^)i/^^du  dv  =  P^  , 
\  U   {-^   +  v^)i^_du  dv  =  P_  . 


5.   Series  Solution  of  the  Collision-free  Boltzmann  Equations. 

Our  procedure  for  solving  the  system  (22),  (23),  (28),  (29), 
with  Initial  conditions  (31)-(56)  Is  to  put 

f  =fi°)  +  .fil)  +e2fp)  +  ... 

+       +  +  + 

in  the  Boltzmann  equations  (28),  (29),  and  equate  the  coefficient 
of  each  power  of  €  on  the  left  with  the  coefficient  of  the  same 
power  of  €  on  the  right  in  the  usual  way.  We  thereby  obtain  a 
sequence  of  equations,  which  we  can  solve  explicitly,  and  we 
obtain  f^°S  f^"^^,  etc.  in  terms  of  B(x),  E(x)  and  their  deriva- 
tives.  A  different,  but  related  scheme,  is  used  by  Chew, 
Goldberger  and  Low  [3],  the  difference  being  in  the  scale 
transformations.  We  then  Insert  our  solution  of  the  Boltzmann 
equation  into  the  Maxwell  equations  (28),  (29)  and  obtain 
differential  equations  --  to  any  order  --  E(x),  B(x).     '' 
Retaining  only  the  terms  of  lowest  significant  order,  we  finally 
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arrive  at  an  equation  of  third  order  for  B(x).  This  equation 
can  be  solved  by  quadratures.   The  solution  yields  the  asymptotic 
shape  of  the  solitary  wave  in  a  warm  plasma. 

To  carry  out  the  iteration  process  for  solving  the  Boltzmann 
equations,  we  can  work  with  the  ion  equation  (22)  only,  since 
the  results  for  the  electron  equation  can  be  obtained  from  the 
result  for  the  ion  equation  by  replacing 

Thus  we  may  consider  the  following: 


(37)         1^  +  ^^BM  §  =  c(i™)  II  -  c2b(x)  31  +  .5a(x)  |i 


v/here 

l/?  l/^-' 

u  =  v;  ^    cos   0,  V  =  w  '^  "sin   9 

'.;ith  the  conditions   for     x  — >  -oo, 

(38)  f(x,u,v)  -^    i^(w) 

(39)  B(x)  -^  0 

(40)  E(x)   ->  0    . 

We   solve    (37)   by    setting 


(^^1)  ^  =  ^o  +  ^^1  ■*"  ^^^2  +  ^^^3  +  ^^^4  + 
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and  solving  the  sequence  of  equations,  corresponding  to  the  suc- 
cessive powers  of  e,  which  we  get  on  Inserting  (^1)  Into  (37). 
We  use  as  Initial  conditions,  to  determine  the  constants  of 
Integration,  the  conditions 
for  X 


(42)      f^-^i',              f-^  ^  0,  f 2  ->  0  ,  etc. 

(4:5)      B(x)  ^^  0,    B'(x)-^0,  B"(x)  -^  0,   etc. 

(44)      E(x)  ->  0,    E'(x)  — $>0,  E"(x)  ->  0,  etc. 

We  shall  have  to  compute  f  ,  f, ,  fp,  f^,  and  f^  In  order  to 
obtain  the  lowest  order  approximation  of  the  field  variables 
B,  etc. 

In  the  lowest  order,  we  have 


O   _  r. 


whence 


^o  =  ^o(^'^) 


I'/here  ^     Is  yet  to  be  determined.  The  function  f6     will  be  deter- 
mined by  a  compatibility  condition,  for  the  solvability  of  the 
next -order  equation  for  f-,.   In  the  next  order,  we  have 

(^5)      ^=  (1+-)^=  (1+u)  J^^^(x,w)  . 
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We  integrate  both  sides  with  respect  to  9,   holding  x,w  fixed, 
from  e  =  0  to  0  =  2ir.      We  get 

(46)         0  =  ^° 


Using  the  initial  condition  (42),  we  get 


(47)  i^Q  =  T^(w);   so  that    f^  =  i^(w). 


Now,  from  (45),  (46)  we  have 

Sf, 


Sf-, 
(48)  1=0 


v/hence 

(49)  fi  =  i2^i(x,w). 

Again,  16-.    will  be  determined  from  the  next  equation  (coefficient 
of  e^). 

This  next  equation  is 

ht  Bf        Sf,     ^f 

Using  (47),  (49) J  this  becomes 

(51)  ^  =  (1+u)  ^  ;2$-l(x,w)  +  BV''(w)[-2v]. 


Integrating  v;ith  respect  to  9,   we  get 


(52)  0=|3^^^(x,w) 


Thus,  by  (42),  we  obtain 
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(55)  )^1  =  0;   so  that   f^  =  0. 

Now  (51)  becomes 

^  =  B^'(w)[-2v] 


whence 

(5^)  fg  =  BT^'(w)[2u]  +  jz52(x,w). 

Again,  to  find  s6^   we  go  to  the  next  equation  (coefficient  of  €  ). 
The  next  equation  is 

Sf,     ^f.        Sfp     Bf,     Sf^ 

(55)  ^  +  B  ^  =  (1+u)  ^  -  B  ^  +  E  ^ 

or,  by  (47),  (53),  (5^), 

Sf^  ^)2^o  P 

(56)  ^=  (1+u)  ^  +  B'V"(w)[2u+2u'^]+E^'(w)[2u]  . 
Integrating  (56)  with  respect  to  9,   we  obtain 

(57)  °  =  B3r  ■*•  B'wi^'(w) 
whence 

i^Jg  =  -  BwV  (w) 
and,  by  (54), 

(58)  fg  =  BTl''(w)[-w  +  2u]. 

Multiplying  (57)  by  (1+u)  and  subtracting  the  result  from  (56), 
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vie   get 

(59)  ^^  =  B'-^'(w)[-w  +  (2-w)u  +  2u  ]  +  Ki/' {vj)[2u]. 

We  can  find  an  integral  of  this  equation:   We  make  use  of  the 
facts 

We  find 

(60)  f^  =  B'T^' (w)[(2-w)v  +  uv]  +  E^'(w)[2v]  +  i2$,(x,w)  . 
Again,  to  find  ^.^   we  must  use  the  next  equation  (coefficient  of 

By  nov;  the  nature  of  the  algorithm  should  be  clear.   From 
this  point  on  the  calculations  rapidly  become  lengthier.   There- 
fore we  shall  not  write  out  all  the  steps,  but  at  this  point 
merely  give  the  results.  We  obtain  the  following: 

(61)  f^  =  if 

(62)  f^  =  0 

(63)  f^   =  BT^'t-^  +  2u] 

(64)  f,  =  B'T^'[(2-w)v  +  uv]  +  Ei^'[2v] 

(65)  f^  =  B"^'[(J  w  -  §  w^)  +  (-2+w)u  +  (-  1  +  i-w)u2-  ^u^] 

+  E't/^'LI  w  -  2u  -  U-]  +  B^T^'[(l+w)-4u] 
+  B^tA"[^  w^  -  2wu  +  2u^]   . 
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4.   Computation  of  the  Moments 

We  now  compute  the  moments  of  f  ,  f_  which  occur  in  the 
Maxwell  equations.  We  have 


(66)      //  •^(w)du  dv  =  TT  /  i^(w)dw  =  1 

o 


and  we  define     P,    6     by 


(67)  2P  =  //  W!^(w)du  dv  =  TT  /     wi^(w)dw 

o 

op  °°      p 

(68)  4dP     =  //  w  i^(w)du  dv  =  tt  /     w  •\fr(w)dw. 


(Note  that    (68)   is  a   fourth  moment.)     By  applying  Schwarz's 
Inequality   (|/F-G|^  <    |/F^|'|/G^|)     it   can  be  shown  that  if  if 
is  positive,   then 

(5  >  1 

and  that  if  if'(w)  is  positive  and  a  monotone  decreasing  function 
of  w,  that 

6   >  4/5. 
If  1^  is  a  Maxwellian  distribution,  then 

6  =  2. 
By  integrating  (61)  to  (65)  we  obtain  the  following: 
(69)      ;/  f^du  dv  =  1 
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(70)  ;/  f3_du  dv  =  0 

(71)  //  fgdu  dv  =  B 

(72)  //  f^du  dv  =  0 

(73)  11  f^du  dv  =  B"[-l  +  I  P]  -  E'  +  2B^. 

Also,  we  obtain  the  following: 

(74)  //  vf^du  dv  =  0 

(75)  //  vf^du  dv  =  0 

(76)  /J  vf^du  dv  =  0 

(77)  //  vf^du  dv  =  B'[-l   +  2P] 

(78)  //  vfi^du  dv   =   0. 

We  shall  also  need  the  moment  of  fr-.   This  can  be 

5 

obtained,  without  actually  computing  f^.,  as  follows:  we 
multiply  the  Boltzmann  equation  (37)  by  u  and  integrate, 
noting  that 

//  u  1^  du  dv  =  ;/  u[u  |i  -  V  |I]  du  dv 

=  //  [^  (u^f)  -  ^  (uvf)]du  dv  +  //  vf  du  dv 

=  //  vf  du  dv. 
We  obtain,  then, 
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(l+eS)  //vf  du  dv  =  €  J^  J/(u+u^)f  du  dv  -  e^E  //f  du  dv 

so  that 

//v  f^du  dv  +  B//vf^du  dv  =  -^  //(u+u^)f2^du  dv  -  E  //fgdu  dv. 

In  this  way  we  obtain 

(79)  //vf  du  dv  =  B"'[l  -  TT  P  +  I  6P^]  +  E"[l  -  |  P]+  3BB' . 
Prom  these  results  we  obtain  the  following: 

(80)  //f^du  dv  =  1  +  eS  +  e^  [b"[-1  +  |  P^]-  E'  +  2B^j  +  O(e^) 


(81)  //f_du  dv  =  1 +€S  +  e^[  a^B"[-l  +  I  P_]+  aE'+  2B^+0(e^) 

(82)  //vf_|_du  dv  =  €^|b'[-1  +  2F^]    -  eJ 

+  €5{b"'[1-  ^  P+  +  I  6+P^]  +  E"[l-  I  P+]+  3BB']+0(e'^) 


(85)      //vf_du  dv  =  €^  j  B'[a-2aP_]-  e) 

+  €5|B'"[-a^+  ^  aP?_-   I  a^6_P^]+  E"[a^-  |  a^Pj 
-  3aBB']  +  0(€'^)  . 

These  moments  are  sufficient  for  our  purpose  in  the  following 
sections.   However,  an  appreciable  amount  of  additional  calculations 
will  be  necessary  in  order  to  obtain  an  explicit  expression  for  the 
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anisotropic  pressure  tensor  corresponding  to  our  distribution 
function.   Then  a  more  direct  comparison  of  our  results  with 
other  two-fluid  theories  will  be  possible. 


5.   Differential  Equation  for  B 

If  we  insert  the  expressions  (80),  (81),  (82),  (8p)  into 
the  Maxwell  equations  (28),  (29)  we  obtain  the  following: 

2 

(84)  (1+a)  ^  E'  =  B"[-(l-a^)  +  |  (P^-a2p_ ]-(l+a)E'+  0{e^) 

c 

(85)  [-(1+a)  +  (l+€^)(2P_j_+  2aP_)]B' 

=  (l+€^)|B'[-(l+a)  +  (2P_^  +  2aP_)]) 

+  e2^B'"[(l+a^)-  ^  (P^  +  a^P_)+   |(6^p2  +  aV)] 

+  E"[(l-a^)-  I  (P_^  -  a^Pj]  +  3(l+a)BB'|  +  0(€^). 

Since  the  lov;est-order  terms  cancel,  (85)  can  be  written  as 
follows : 

(86)  B"' [(1+a^)-  ^  [P^  +  a^Pj  +  |  (C^P^  +  a^d_P?)] 

+  E"[(l-a^)-  I  (P^-a2p_)]  +  :5(l+a)BB'-  (l+a)B'  =  0(€^). 

Now  we  eliminate  E  from  (86)  by  use  of  (84).   In  this 
way  we  find  a  differential  equation  for  B(x)  of  the  following 
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2 

form  (neglecting  terms  which  are  0(e  )). 

(87)      KB'"  +  5BB'  -  B'  =0. 

The  coefficient  R  is  quite  lengthy  but  can  be  simplified  some- 
what  if  the  term  containing  A  /c   in  (84)  is  neglected.   Neglecting 
this  term  is  equivalent  to  ass\jm.ing  charge  neutrality,  instead  of 
using  Polsson's  equation.   With  this  simplification  the  coefficient 
R  becomes 


(88)      R  =  a  +  i^-^^"^)    (p  -  a^p  ) 

4(l+a)    ^ 


2„  n2 


+  i-^   6(l+a)(e,P^  +  a^6  p2)-  9(P  -  a^P  ) 

4(l+a)'^  I        "^"^      '"      ^  ~      ^ 

It  is  convenient  to  rewrite  (88)  in  terms  of  p  ,  p_,  which 
are  defined  as  the  ratios  of  the  pressures  of  the  ion  gas  and 
electron  gas,  respectively,  to  the  magnetic  pressure.   That  is. 


P.  P 


+   B^/2m.    '       -   B^/2|a 


We  have  also 


(l+a)6  (l+a)p_ 

(90)      P^  =  — 5— ^t ,   aP_  = 


■+   2(1+€^)(1+P^+P_)  '  ~        2(1+€^)(1+P^+P_) 

Inserting  these  expressions  into  (88),  we  obtain 

(l-12a)(6  -aP  ) 

(91)      R  =  a  + 

8(H-p_^+p_) 

+  l^|6(l+a)(Q^+  ^l  +  a6_^l)   -   90^.-aP_)^j   . 
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We  may  assume 

a  «  1 

since  the  mass  of  an  electron  Is  much  less  than  the  mass  of  any 
positive  ion.   If  we  also  assume  that  the  ion  and  electron 
temperatures  are  equal,  so  that  the  pressures  are  equal,  we 
have 

For  small  p,  which  is  an  interesting  case,  (91)  reduces  to 
(92)  R  =  a  +  ^  . 

Comparing  this  with  the  expression  we  get  if  we  do  not  neglect 

2  /  o 
A  /c  ,  we  see  that  a  necessary  condition  for  the  validity  of 

the  charge -neutral  assumption  is 
(95)  4«  ("  +  &)• 

ill 
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6.   width  of  the  Solitary  Wave 

The  differential  equation  (87)  can  be  solved  by  quadratures. 
The  solution  satisfying  the  condition 

for  X  — >  -00  ,    B  — >  0 
Is 

(94)  B(x)  = 


cosh^^  ^ 


2R^ 

(apart  from  a  constant  which  may  be  added  to  x;  this,  however, 
may  be  removed  by  a  shift  of  coordinates). 

Plotted  against  x,  this  is  a  symmetric  curve,  rising  from 
B  =  0  at  X  =  -00  to  its  maximum  B  =  1  at  x  =  0,  and  then  again 
falling  to  B  =  0  at  X  =  +00  .  And  B  takes  on  half  its  maximum 
value  at  the  two  points  where 


--  i^) '  ^^'^ 


that  is,  at  the  two  points 

X  =  t  2R-^/^x  0.881   . 

Now  by  (14),  (15),  (16)  defining  the  dimensionless  variables  x, 
we  obtain  from  this  the  following  expression  for  the  width  L  of 
the  solitary  wave : 

-.  -^    /     \       rfi   *  v/m  m  "  U 


W- 


o 
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which  should  be  accurate  for  small  (positive)  values  of  (m  -  l), 
and  small  p,  if  the  condition  (95)  for  the  validity  of  the 
charge -neutral  assumption  obtains.   Of  course,  under  these 
assumptions  U   is  nearly  the  Alfven  speed,  so  we  have  the 
following  alternative  for  L: 


m.   ■   -  /  m  ^ 


L^-1^52_  J-L  ^  1  (_+)     1 

^^-^— n  /    IF  VmJP   e  /ixn^ 

The  analogous  formula  for  the  cold  plasma  --  which  can  also 
be  obtained  directly  (our  earlier  work)  --  is  found  by  setting 
p  =  0. 

Thus  we  see  that  the  effect  of  plasma  temperature  on  the 
Adlam-Allen  solitary  wave  is  to  multiply  the  width  of  the 
wave  by  the  factor 


/ 


1  +  TO  (^)p 


for  small  g. 
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